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Abstract 

Using 1.0 fb^^ of pp collision data collected at a centre-of-mass energy of ^/s = 7 TeV 
with the LHCb detector, measurements of the polarization amplitudes, strong phase 
difference and triple product asymmetries in the — )• (^(/> decay mode are presented. 
The measured values are 

l^ol^ = 0.365 ± 0.022 (stat) ± 0.012 (syst) , 
= 0.291 ± 0.024 (stat) ± 0.010 (syst) , 
cos(5||) =-0.844 ± 0.068 (stat) ± 0.029 (syst) , 
Au =-0.055 ± 0.036 (stat) ± 0.018 (syst) , 
Av = 0.010 ± 0.036 (stat) ± 0.018 (syst) . 
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1 Introduction 



In the Standard Model, the flavour- changing neutral current decay — t- (jxf) proceeds 
via a 6 ^ sss penguin process. Studies of the polarization amplitudes and triple product 
asymmetries in this decay provide powerful tests for the presence of contributions from 
processes beyond the Standard Model flfls]. 



The (j)(j) decay is a pseudoscalar to vector- vector transition. As a result, there are 

three possible spin configurations of the vector meson pair allowed by angular momentum 
conservation. These manifest themselves as three helicity states, with amplitudes denoted 
H^i,H_i and Hq. It is convenient to define linear polarization amplitudes, which are 
related to the helicity amplitudes through the following transformations 



All 



Ho 
K 



+1 



V2 



(1) 



The (f)(f) final state can be a mixture of CP-even and CP-odd eigenstates. The lon- 
gitudinal {Aq) and parallel (Ay) components are CP-even and the perpendicular com- 
ponent {A±) is CP-odd. From the V-A structure of the weak interaction, the longitu- 
dinal component, fi = \Aq\^/{\Aq\^ + -f- is expected to be dominant [6[]8]. 



However, roughly equal longitudinal and transverse components are found in measure- 
ments of B+ (pK*+, 5° ^ (pK*^, B+ p°K*+ and S° ^ p°K*° decays at the 



B-factories [9 14 . To explain this, large contributions from either penguin annihilation 
effects [15] or final state interactions 16 have been proposed. Recent calculations where 



phenomenological parameters are adjusted to account for the data allow fi in the range 
0.4 — 0.7 ^[t]. Another pseudoscalar to vector- vector penguin decay is B^ — )■ K*^K*^. A 
recent measurement by the LHCb Collaboration in this decay mode has found a value of 
h = 0.31 ±0.12 ±0.04 (17. 

The time-dependent differential decay rate for the 5° — )■ (p(f) mode can be written as 



d cos 9id cos 92d^dt 



(2) 



i=l 



where the helicity angles n = (9i, 62, $) are defined in Fig. [Tj The angular functions 
are flsl 



/3(^1,^2,$ 
/4(^1,^2,$ 



4 cos^ 9i cos^ 6^2 , 
sin 9 1 sin^ ^2(1± cos 2$) 
sin^ 9i sin^ ^2(1 - cos 2$) 
-2 sin^ 9i sin^ ^2 sin 2$ , 
72 sin 2^1 sin 2^2 cos $ , 
sin 29i sin 2^2 sin $ . 



(3) 



1 




Figure 1: Decay angles for the — )■ 00 decay, where the K'^ momentum in the 0i^2 rest 
frame, and the parent 0i,2 momentum in the rest frame of the meson span the two 
meson decay planes, 61^2 is the angle between the track momentum in the 0i^2 nieson 
rest frame and the parent 0i^2 momentum in the B^ rest frame, $ is the angle between 
the two meson decay planes and ni 2 is the unit vector normal to the decay plane of the 
01,2 meson. 



The time-dependent functions Ki{t) are given by 19 
1 



Ki{t) = ^A^[(l + cos0s)e-^L* + (l-cos0,)e~^H*±2e-^=*sin(Am,t)sin. 
1 



K2{t) = + cos0s)e-^Lt ^ (X _cos0,)e"^H*±2e"^=*sin(Amst)sin. 

1 



K^it) = ^Al [(1 - cos 0s)e"'^Lt + (1 + cos 0s)e"'^«* T 2e"^=* sin(Amst) sin 0^] , 
K4{t) = |A||||A_L|[±e^^°*{sin5i cos(Amst) — cos(5i sin(Amst) cos0s} 
-^(e^^H* - e~^i^*) cos5i sin 0s] , 

K,{t) = ^|Aopj||cos(52-5i) 

[(1 + cos 0s)e~^L* + (1 - cos 0s)e"^Ht _^ 26"^-* sin(Amst) sin 0^] , 
Kg^t) = |y4o| |v4_l| [±e~^°*{sin (^2 cos(Amst) — cos ^2 sin(Amst) cos0s} 

(e-^«* - e-^i^*) cos 62 sin 0^] , (4) 



1 

2' 



where the upper of the ± or =p signs refers to the B^ meson and the lower refers to 
a 5° meson. Here, Fl and Fh are the decay widths of the light and heavy B^ mass 
eigenstatesj^ Ams is the B^ oscillation frequency. Si = aTg{A±/ A\\) and S2 = 8j:g{A±/ Aq) 
are CP-conserving strong phases and 0^ is the weak CP-violating phase. It is assumed 



^Units are adopted such that h = 1. 



2 



that the weak phases of the three polarization amphtudes are equaL The quantities Fh 
and Fl correspond to the observables AF^ = Fl — Fh and F^ = (Fl + Fh)/2. In the 
Standard Model, the value of (ps for this mode is expected to be very close to zero due 
to a cancellation between the phases arising from mixing and decay 20 ,^1A calculation 
based on QCD factorization provides an upper limit of 0.02 rad for 0^ [6,21 . This is 
different to the situation in the — )• J/ipcp decay, where the Standard Model predicts 



(/),(J/^0) = -2aTg{-VtsVtl/VcsV;^) = -0.036 ± 0.002 rad 22 . The magnitude of both 



weak phase differences can be enhanced in the presence of new physics in 5° mixing. 



where recent results from LHCb have placed stringent constraints 23 . For the — )■ (fxf) 
decay, new particles could also contribute in 6 — )■ s penguin loops. 

To measure the polarization amplitudes, a time-integrated untagged analysis is per- 
formed, assuming that an equal number of and B^ mesons are produced and that the 
CP-violating phase is zero as predicted in the Standard Modelj^ In this case, the functions 
Ki(t) integrate to 





= l^olVrL, 


K2 


= l^iilVrL, 




= l^±r/rH, 




= 0, 


K, 


= v4|| cos 




= 0, 



(5) 

where the strong phase difference is defined by d\\ = 62 — Si = arg(y4||/y4o) and the time 
integration assumes uniform time acceptance. 

In addition, a search for physics beyond the Standard Model is performed by studying 
the triple product asymmetries [l}]3] in the B'^ — )■ (fxj) decay. Non-zero values of these 
quantities can be either due to T-violation or final-state interactions. Assuming CPT 
conservation, the former case implies that CP is violated. Experimentally, the extrac- 
tion of the triple product asymmetries is straightforward and provides a measure of CP 
violation that does not require flavour tagging or a time-dependent analysis. 

There are two observable triple products denoted U = sin(2$)/2 and V = ±sin($), 
where the positive sign is taken if the T-even quantity cos 9i cos 6*2 > and the negative 
sign otherwise. These variables correspond to the T-odd triple products 

sin $ = {fii X 722) ■ pi , 
sin(2$)/2 = (hi ■ n2){ni X 712) ■ pi , (6) 

where fii {i = 1,2) is a unit vector perpendicular to the (pi decay plane and pi is a unit 
vector in the direction of the 0i momentum in the B^ rest frame. The triple products, U 

■^The convention used in this Letter is that the symbol (/j^ refers solely to the weak phase difference 
measured in the i3° — >■ decay. 

"^In the case of non-zero 0s deviations from these formulas are suppressed by a factor of AFs/r^ and 
hence only small variations would be observed on the fitted parameters. 
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and V, are proportional to the /4 and /g angular functions which, for 0^ = 0, vanish in 
the untagged decay rate for any value of t. The and /e angular functions would not 
vanish in the presence of new physics processes that cause the polarization amplitudes 
to have different weak phases [l]. Therefore, a measurement of significant asymmetries 
would be an unambiguous signal for the effects of new physics (l]|3]. 
The asymmetry, Au, is defined as 

where (A^-) is the number of events with f/ > (f/ < 0). Similarly Ay is defined as 

where M+ (M_) is the number of events with > < 0). The triple product asym- 
metries, Au and Ay are proportional to the interference terms Xm{A±A*^) and Xm{A±AQ) 
in the decay rate. 

The Bg — > (f)(j) decay mode was first observed by the CDF Collaboration [24|. More re- 
cently, CDF has reported measurements of the polarization amplitudes and triple product 



asymmetries in this mode based on a sample of 295 events 25 . In this Letter, measure- 
ments of the polarization amplitudes, |AoP and |A±P, the strong phase difference, S\\, 
and the triple product asymmetries, Au and Ay, are presented. The dataset consists of 
801 ± 29 candidates collected in 1.0 fb~^ of pp collisions at the LHC. The Monte Carlo 
(MC) simulation samples used are based on the Pythia 6.4 generator [26] configured 



with the parameters detailed in Ref. 27 . The EvtGen 28 and Geant4 m)\ packages 



are used to generate hadron decays and simulate interactions in the detector, respectively. 



2 Detector description 

The LHCb detector ^O] is a single-arm forward spectrometer covering the pseudorapid- 
ity range 2 < r/ < 5, designed for the study of particles containing 6 or c quarks. The 
detector includes a high precision tracking system consisting of a silicon-strip vertex de- 
tector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift-tubes placed downstream. The combined tracking 
system has a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at 
100 GeV/ c, and an impact parameter resolution of 20 [im for tracks with high transverse 
momentum. Charged hadrons are identified using two ring-imaging Cherenkov detectors. 
Photon, electron and hadron candidates are identified by a calorimeter system consist- 
ing of scintillating-pad and pre-shower detectors, an electromagnetic calorimeter and a 
hadronic calorimeter. Muons are identified by a muon system composed of alternating 
layers of iron and detector stations. The trigger consists of a hardware stage, based on 
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information from the calorimeter and muon systems, followed by a software stage which 
applies a full event reconstruction. 

The software trigger used in this analysis requires a two-, three- or four-track sec- 
ondary vertex with a high sum of the transverse momentum, p^, of the tracks, significant 
displacement from the primary interaction, and at least one track with > 1.7GeV/c; 
impact parameter with respect to the primary interaction greater than 16; and a track 
fit x^/ndf < 2 where ndf is the number of degrees of freedom in the track fit. A multivari- 
ate algorithm is used for the identification of the secondary vertices [31j. The — (pep 
candidates are selected with high efficiency either by identifying events containing a (p 
meson or using topological information to select hadronic b decays. Events passing the 
software trigger are stored for subsequent offiine processing. 



3 Event selection 

The Bg — 7- (j)(f) channel is reconstructed using events where both mesons decay into a 
K^K~ pair. The 5^ — t- cfxp selection criteria were optimized using a data-driven approach 
based on the gVlot technique employing the four-kaon mass as the unfolding variable |32| 
to separate signal (S) and background (B) with the aim of maximizing S/ y/S + B. The 
resulting cuts are summarized in Table [Tj Good quality track reconstruction is ensured 
by a cut on the transverse momentum (px) of the daughter particles and a cut on the 
X^/ndf of the track fit. 

Combinatorial background is reduced by cuts on the minimum impact parameter sig- 
nificance of the tracks with respect to all reconstructed pp interaction vertices and also by 
imposing a requirement on the vertex separation oi the B^ candidate. Well-identified (p 
meson candidates are selected by requiring that the two particles involved are identified as 
kaons by the ring-imaging Cherenkov detectors using a cut on the difference in the global 
hkelihood between the kaon and pion hypotheses (A In Ck-k > 0) and by requiring that 
the reconstructed mass of each K^K~ pair is within 12 MeV/c^ of the nominal mass of 



the (p meson |33 . Further signal purity is achieved by cuts on the transverse momentum 
of the (p candidates. 

Figure [2] shows the four-kaon invariant mass distribution for selected events. To de- 
termine the signal yield an unbinned maximum likelihood fit is performed. The -B° — )■ cpcp 
signal component is modelled by two Gaussian functions with a common mean. The res- 
olution of the first Gaussian is measured from data to be 13.9 ± 0.6 MeV/c^. The relative 
fraction and resolution of the second Gaussian are fixed to 0.785 and 29.5 MeV/c^ respec- 
tively, where values have been obtained from simulation. Combinatorial background is 
modelled using an exponential function. Background from B'^ — (pK*^ and 5° — )■ K*^K*^ 
decays is found to be negligible both in simulation and data driven studies. Fitting the 
probability density function (PDF) described above to the data, a signal yield of 801 ± 29 
events is found. 

In addition to the dominant P-wave (p — )■ K^K^ component described in Section [l] 
other contributions, either from /q — > K^K^ or non-resonant K~^K~, are possible. The 
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Table 1: Selection criteria for the B'^ — )■ 00 decay. The abbreviation IP stands for impact 
parameter and pt*^^ and pt'^'^ refer to the transverse momentum of the two candidates. 



Variable 


Value 


Track x^/ndf 


< 


5 


Track 


> 


500 MeV/c 


Track IP 


> 


21 


A In Ckvt 


> 







< 


12 MeV/c2 




> 


900 MeV/c 
2 GeVVc^ 




> 


vertex x^/ndf 
B^ vertex x^/ndf 


< 


24 


< 


7.5 


5° vertex separation 


> 


270 


5° IP 


< 


15 



size of these contributions, neglecting interference effects, is studied by relaxing the 
mass cut to be within 25 MeV/c^ of the nominal valu^ and using the gVlot technique in 
conjunction with the mass to subtract the combinatorial background. 

The resulting mass distribution is shown in Fig. [3j A fit of a relativistic P-wave 
Breit-Wigner function together with a two body phase space component to model the 
S-wave contribution is superimposed. In a ±25 MeV/c^ mass window, the size of the 
S-wave component is found to be (1.3 ± 1.2)%. Since the S-wave yield is consistent with 
zero, it will be neglected in the following section. A systematic uncertainty arising from 
this assumption will be assigned. 

4 Results 

The polarization amplitudes (|v4op, |^±P, are determined by performing an un- 

binned maximum likelihood fit to the reconstructed mass and helicity angle distributions. 
For each event, the meson used to define 9i is chosen at random. Both the signal and 
background PDFs are the products of a mass component described in Section [3] together 
with an angular component. The angular component of the signal is given by Eq. [3] mul- 
tiplied by the angular acceptance of the detector. The acceptance is determined using 
the simulation and is calculated separately according to trigger type, i.e. whether the 
event was triggered by the signal candidate or other particles in the event. In total the fit 
for the polarization amplitudes has eight free parameters: the signal angular parameters 
|^op5 l^-Lp and cos((5|j) defined in Section [l| the fractions of signal for each trigger type, 

^This is a larger window than the ±12 MeV/c^ window used in the polarization amplitude and strong 
phase difference measurements. 
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Figure 2: Invariant K~^K^K^K^ mass distribution for selected — )■ 00 candidates. 
A fit of a double Gaussian signal component together with an exponential background 
(dotted line) is superimposed. 



the resolution of the core Gaussian, the mass and the slope of the mass background. 

„^ ...... ..... , . 12 , 1/1. 12 , I4|j| 



The sum of squared amplitudes is constrained such that l^ol + I^J-I + l^lll = 1- The 
angular distributions for the background have been studied using the mass sidebands in 
the data, where mass sidebands are defined to be between 60 and 300 MeV/c^ either side 
of the nominal 5° mass fSS]. With the current sample size these distributions are consis- 
tent with being flat in (cos 6'i, cos 6*2, $). Therefore, a uniform angular PDF is assumed 
and more complicated shapes are considered as part of the systematic studies. The values 
of r, = 0.657 ± 0.009 ± 0.008 ps-i and AF, = 0.123 ± 0.029 ± 0.011 ps^^ together with 
their correlation coefficient of —0.3 quoted in 23 are used as a Gaussian constraint. The 
validity of the fit model has been extensively tested using simulated data samples. The 
results are given in Table |2] and the angular projections are shown in Fig. |4j 

Several sources of systematic uncertainty on the determination of the polarization am- 
plitudes are considered and summarized in Table |3| With the present size of the dataset, 
the S-wave component is consistent with zero. From the studies described in Section [3] 
and fits to the data including the S-wave terms in the PDF ^4], we consider a maximum 
S-wave component of 2%. Simulation studies have been performed to investigate the effect 
of neglecting an S-wave component of this size. As discussed in Section [T| the integration 
that leads to Eq. [5] assumes uniform time acceptance. This is not the case due to lifetime 
biasing cuts in the trigger and offline selections. The functional form of the decay time 
acceptance is obtained through the use of Monte Carlo events. The difference between us- 
ing this functional form in simulation studies and using uniform time acceptance is taken 
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1000 1010 1020 1030 1040, 

(MeV/c^) 

Figure 3: Invariant mass distribution of K^K^ pairs for the 1?° — )■ 00 data without a 
mass cut. The background has been removed using the ^P/ot technique in conjunction 
with the K^K~ invariant mass. There are two entries per 5° candidate. The sohd hne 
shows the resuh of the fit model described in the text. The fitted S-wave component is 
shown by the dotted hne. 

Table 2: Measured polarization amplitudes and strong phase difference. The uncertainties 
are statistical only. The sum of the squared amplitudes is constrained to unity. The 
correlation coefficient between \Aq\^ and \A±_\^ is —0.47. 



Parameter 


Measurement 




0.365±0.022 




0.291±0.024 


|yl|ll' = l-(|Aop + |Axp) 


0.344±0.024 


C0S((5||) 


-0.844±0.068 



as a systematic uncertainty. The uncertainty on the angular acceptance for the signal is 
propagated to the observables also using Monte Carlo studies. The analysis was repeated 
with an alternative background angular distribution, taken from a coarsely binned his- 
togram in (cos 6'i, cos 92-,^) of the mass sidebands, and the difference taken as a systematic 
uncertainty. An additional uncertainty arises from angular acceptance dependencies on 
trigger type. This dependency is corrected for using Monte Carlo events, with half of 
the effect on fitted parameters assigned as systematic uncertainties. The total systematic 
uncertainty is obtained from the sum in quadrature of the individual uncertainties. 
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Figure 4: Angular distributions for (a) $, (b) cos 6*1 and (c) cos 6*2 of — )■ 00 events 
with the fit projections for signal and background superimposed for the total fitted PDF 
(solid line) and background component (dotted line). 

Table 3: Systematic uncertainties on the measured polarization amplitudes and the strong 
phase difference. 



Source 


1 A |2 
l^ol 






cos 5|| 


S-wave component 


0.007 


0.005 


0.012 


0.001 


Decay time acceptance 


0.006 


0.006 


0.002 


0.007 


Angular acceptance 


0.007 


0.006 


0.006 


0.028 


Trigger category 


0.003 


0.002 


0.001 


0.004 


Background model 


0.001 




0.001 


0.003 


Total 


0.012 


0.010 


0.014 


0.029 



The distributions of the U and V triple product observables are shown in Fig. |5] for 
the mass range 5286.6 < M{K+K-K+K-) < 5446.6 MeV/c^. To determine the triple 
product asymmetries, the dataset is partitioned according to whether U (V) is less than 
or greater than zero. Simultaneous fits are performed to the mass distributions for each 
of the two partitions corresponding to each observable individually. In these fits, the 
mean and resolution of the Gaussian signal component together with the slope of the 
exponential background component are common parameters. The asymmetries are left as 
free parameters and are fitted for directly in the simultaneous fit. The measured values 
are 

Au = -0.055 ±0.036 , 
Av = 0.010 ±0.036 . 

Systematic uncertainties due to the residual effect of the decay time, geometrical accep- 
tance and the signal and background fit models have been evaluated and are summarized 
in Table |4} The effect of the decay time acceptance has been found using the same method 
as for the polarization amplitudes. The impact of angular acceptance on the measured 
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Figure 5: Distributions of the U (left) and V (right) observables for the — )■ (jxp data in 
the mass range 5286.6 < M{K+K~K+K-) < 5446.6 MeV/c^ The distribution for the 
background is taken from the mass sidebands, normahzed to the same mass range and is 
shown by the sohd histogram. 



values has been obtained from simplified simulation studies. The total systematic uncer- 
tainty is conservatively estimated by choosing the larger of the two individual systematic 
uncertainties on Au and Ay. The contributions are combined in quadrature to determine 
the total systematic error. Various cross-checks of the stability of the result have been 
performed. For example, dividing the data according to how the event was triggered or 
by magnet polarity. No significant bias is observed in these studies. 

Table 4: Systematic uncertainties on the triple product asymmetries Ajj and Ay- The 
total uncertainty is the quadratic sum of the larger of the two components. 



Source 


Au 


Av 


Final uncertainty 


Angular acceptance 


0.009 


0.006 


0.009 


Decay time acceptance 


0.006 


0.014 


0.014 


Fit model 


0.004 


0.005 


0.005 


Total 


0.018 



5 Summary 

The polarization amplitudes and strong phase difference in the (f)(j) decay mode are 

measured to be 

0.365 ± 0.022 (stat) ± 0.012 (syst) , 
0.291 ± 0.024 (stat) ± 0.010 (syst) , 
0.344 ±0.024 (stat) ± 0.014 (syst) , 
0.844 ± 0.068 (stat) ± 0.029 (syst) , 



\Ao? = 
\A±\' = 

cos(5||) = - 
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where the sum of the squared amphtudes is constrained to be unity. These values agree 
well with the CDF measurements 2 

It is interesting to 



Measurements in other B — )■ VV penguin tran- 
sitions at the B factories generally give higher values of /l (9 14 

note that the value of fi found in the B^ — t- 00 channel is almost equal to that in the 

The results are in agreement with QCD factorization predic- 
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S° K*^K*'^ decay 

tions (6|[7], but disfavour the pQCD estimate given in [s]. 

The triple product asymmetries in this mode are measured to be 

Au = -0.055 ± 0.036 (stat) ±0.018 (syst) , 
Av = 0.010 ±0.036 (stat) ±0.018 (syst) . 

Both values are in good agreement with those reported by the CDF Collaboration 
and consistent with the hypothesis of CP conservation. 
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